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ABSTRACT 

We use the suite of simulations presented in iPoole et al.l (|2006[ ) to examine global 
X-ray and Sunyaev-Zel'dovich (SZ) observables for systems of merging relaxed X-ray 
clusters. The time evolution of our merging systems' X-ray luminosities, tempera- 
tures, total mass measures, SZ central Compton parameters and integrated SZ fluxes 
are presented and the resulting impact on their scaling relations examined. In all 
cases, and for all parameters, we observe a common time evolution: two rapid tran- 
sient increases during first and second pericentric passage, with interceding values 
near or below their initial levels. T his is in good qualitativ e agreement with previous 
idealized merger simulations {e.g. Ri cker fc Sarazinll2001l ), although we find several 
important differences related to the inclusion of radiative cooling in our simulations. 
These trends translate into a generic evolution in the scaling-relation planes as well: 
a rapid transient roughly along the mass scaling relations, a subsequent slow drift 
across the scatter until virialization, followed by a slow evolution along and up the 
mass scaling relations as cooling recovers in the cluster cores. However, this drift is not 
sufficient to account for the observed scatter in the scaling relations. We also study 
the effects of mergers on several theoretical temperature measures of the intracluster 
m edium: emission weight ed measures (T ew ), the spectroscopic- like measure proposed 
bv lMazzotta et al] (|2004L T s /) and plasma model fits to the integrated spectrum of the 
system (T spec ). We find that T s i tracks T spec for the entire duration of our mergers, 
illustrating that it remains a good tool for observational comparison even for highly 
disturbed systems. Furthermore, the transient temperature increases produced during 
first and second pericentric passage are 15 — 40% larger for T ew than for T s i or T spec . 
Th is suggests tha t the e ffects of transient temperature increases on erg and CIm derived 
by iRandall et all (2 002) are over e stimated. Lastly, we examine the X-ray SZ proxy 
proposed by iKravtsov et al] (|2006l ) and find that the tight mass scaling relation they 
predict remains secure through the entire duration of a merger event, independent of 
projection effects. 

Key words: cosmology: theory - galaxies: clusters: general - intergalactic medium 
- X-rays: general 



1 INTRODUCTION 

It has been known for some time that clusters obey simple 
and relatively well defined power-law scalings between many 
of their globally integrated or averaged properties (e.g. tem- 



perature, X-ray luminosity and ma ss). Due to their proven 
utility in cosmological studies (see IVoftll2005l . for an excel- 
lent review), their u se as constraints on theoretical models 
of galaxy formation (|Babul et. al. Il2002f ) and the significant 
discrepancy between o bservations and simple early th eoreti- 
cal expectations (e.g. Kaiser 1986; Markcvitch 1998), these 
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relations have received a great deal of scrutiny in recent 
years. 

Significant progress has recently been made to- 
wards reconciling theoretical expectations for the slope 
and n ormalization of the s e relations with observations 
{e.g. iBoreani et all 12004 IVoit et all 12002ft . As a result, 
attention is now turning towards accounting for their ob- 
served scatte r. One significant early s tep towards this goal 
was made by iMcCarthv et. al. ] (|2004 M04 hereafter) who 
illustrated that the position of systems relative to the mean 
observed L x — T x and Mt — L x relations correlates with 
the system's central entropy, and hence the morphology of 
a system's core. Specifically, they found that systems with 
compact cool cores (CCCs; classically identified as "cooling 
flow" systems) and low central entropies tend to lie on the 
high-luminosity/low-temperature sides of observed distribu- 
tions while systems lacking compact cool cores (NCCs; clas- 
sically identified as "non-cooling flow" systems) and possess- 
ing elevated central ent ropies lie on th e low-l uminosity/low- 
niass sides. Recently, lO'Hara et all {2006) have demon- 
strated that these two populations are indistinguishable in 
these planes when their cores are removed during analy- 
sis. This suggests that variations in core properties are pri- 
marily responsible for the scatter in o bserved scalin g rela- 



tion s, supporting prev ious findings of iFabian et all {l994T ) 
and lMarkevitchl (| 1998ft . 

These results are corroborated by IBalogh et ail (l2006ft 
who find that variations in dark matter structure and un- 
certainties in cosmological models can not account for the 
scatter in the L x —T x or Mt — L x scaling relations. They sug- 
gest that variations in the minimum level of core entropy due 
to differences in heating and/or cooling efficiencies of core 
material are responsible. Little speculation is made regard- 
ing the physical origins of such efficiency variations but they 
theorize that mergers are not likely the source, citing qual- 
itative observations made by authors of previous theoreti- 
cal merger studies that mergers generally d rive systems to 
evolv e parallel to observed scaling relations dRowlev et. al. 
2004). This lies in stark contrast to the work of lSmith et al 
(2005) however, who use strong lensing measurements and 
high resolution Chandra observations to illustrate a correla- 
tion between the degree to which local cluster cores are dis- 
turbed and their variance from the mean Mt — L x , Mt — Tx, 
and L x — T x scaling relations. Hence, several unresolved is- 
sues regarding the role of mergers in shaping the scatter in 
observed cluster scaling relations persist. 

Scaling relations involving the Sunyuev-Zeldovich (SZ) 
properties of clusters, a lthough previously con sidered 
{McCarthy et al.l l2003al lbl; ICavaliere fc Mencil 1200 lft , have 
received significantly less attention than those involving 
mass, X-ray luminosity and temperature. With the impend- 
ing availability of large and statistically significant datasets 
of SZ observations, a careful examination of the effects of 
mergers on SZ scaling relations is warranted. In this paper 
we will use a set of nine ideali zed two-body clus ter merger 
simulations (first presented in IPoole et al. 120061 . Poole06a 
henceforth) to examine the effects of mergers on cluster scal- 
ing relations (refer to Poole06a for a review of pioneerin, 
work involving such simulations and lRicker fc Sarazinll200 
RS01 hereafter, for a particularly good first examination of 
some of the issues revisited in this work). Our simulations 
have been constructed to probe a representative range of im- 



pact parameters and mass ratios of mergers between CCC 
systems with a primary (most massive system) of 10 15 M Q . 
A detailed description of these simulations can be found in 
Poole06a. 

In Section [2] we will briefly summarize several aspects 
of our method for constructing cluster mergers presented 
in Poole06a which are relevant to this paper. In Section [3] 
we present the temporal evolution and scaling relation of our 
systems' X-ray luminosities {L x ) and temperatures {T x ), ex- 
amining the effects of various temperature measures and of 
core excision on the behaviour of these quantities. In Sec- 
tion 2] we examine the effectiveness of hydrostatic estimates 
of total mass {Mt), its isothermal /3-model implementation 
and the M t — T x and M t — L x scaling relations which 
are obtained from them. In Section \5\ we study the evolu- 
tion of our systems' SZ properties (both central Compton 
parameters and integrated SZ fluxes), examining the effec- 
tiveness of the /3-model approximation often implemented 
by observers as well as the most relevant scaling relations 
which can be constructed from them. In Section [6] we ex- 
amine the e ffects of mergers on the X-ray SZ proxy recently 
proposed bv lKravtsov et all {2006). In Section [7] we discuss 
several issues raised by our work and finally, we summarize 
the most notable results of our study in Section [S] 

In all cases our assumed cosmology will be 
(n M ,«A)=(0.3,0.7) with H =75km/s/Mpc. 



2 SIMULATIONS 

We run our simu l ations with GASOLINE 
{Wadslev. Stadel. fc Quinnl |2004| ). a versatile parallel 
SPH tree-code with multi-stepping. We include the effects 
of radiative cooling in our simulations and feedback from 
star formation but the effects of jets from active galactic 
nuclei (AGN) are omitted. 

The basis of our study is a set of 9 idealized 2-body 
cluster merger simulations. Three mass ratios are examined 
(1:1, 3:1 and 10:1) with the most massive system set to have 
a mass of 10 15 Mq in every case. For each mass ratio, three 
impact parameters are examined (one head-on and two off- 
axis cases). We express these impact parameters with the 
quantity v t /V c , where vt is the relative tangential velocity of 
the interacting systems when the core of the smaller system 
crosses the virial radius of the larger system and V c is the 
larger system's circular velocity at its virial radius. 

In Poole06a we present a detailed account of how our 
isolated clusters and their orbits are initialized as well as the 
numerical methods and code parameters of our simulations. 
In this section we summarize key aspects of our approach 
pertinent to our present analysis and direct readers looking 
for additional details to Poole06a. 



2.1 Initial conditions 

We have chosen to focus our study on systems initialized to 
posses compact cores (r c ~ 50 kpc) . We have chosen this ini- 
tial configuration because it is the natural convergent state 
for a relaxed cluster in the absence of feedback and rep- 
resentative of the s tructure of a lar ge fraction of clusters 
{Peres et. al. Ill998l : lEdge et all 1 1992ft . In the context of hi- 
erarchical clustering, such systems are well understood: they 
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Table 1. Percentage difference of the global properties of our systems from the values required to preserve accepted mass scaling 
relations at the end of our simulations and at t re i ax (in parenthesis). For X-ray obscrvablcs, observed mass scaling relations are 
used while for SZ observables (y and Su), mass scaling relati ons computed from our fiducial theoretical models are assumed (see 
Sectionl3.3ll. For temperatur e we use the mass scaling relation of iHorner 1 11200 lh and for luminosity we use the mass scaling relation of 
iReiprich fc Bohrineerl d2002h Corrected results are for analysis conducted with the central O.I-R200 excised and uncorrected values are 
performed without excising the central regions. 



M P \M S vt/v c L x (uncorrected) L x (corrected) T x (uncorrected) T x (corrected) y ^v(-f?250o) 



1:1 


0.00 


21.1 


(-13.1) 


-13.0 


(-2.0) 


10.3 


(5.5) 


2.9 


(-1.5) 


28.8 


(-3.3) 


13.4 


(14.7) 


1:1 


0.15 


-25.6 


(-38.3) 


-13.3 


(-0.6) 


8.2 


(13.1) 


-1.3 


(-0.4) 


-8.7 


(-27.6) 


10.0 


(20.8) 


1:1 


0.40 


-49.6 


(-52.9) 


-27.6 


(-26.1) 


7.5 


(9.2) 


-3.1 


(-4.9) 


-24.2 


(-39.3) 


3.9 


(6.7) 


3:1 


0.00 


-13.3 


(-27.1) 


3.3 


(14.6) 


10.9 


(14.7) 


2.1 


(2.2) 


-17.2 


(-37.7) 


12.9 


(24.7) 


3:1 


0.15 


-38.6 


(-40.4) 


-3.2 


(5.7) 


14.6 


(13.7) 


0.9 


(0.3) 


-29.8 


(-47.6) 


11.8 


(18.2) 


3:1 


0.40 


-31.7 


(-35.2) 


-24.9 


(-16.0) 


5.3 


(9.6) 


-1.7 


(0.0) 


-15.0 


(-33.6) 


-6.6 


(1.6) 


10:1 


0.00 


5.7 


(12.3) 


-10.7 


(19.1) 


3.7 


(13.1) 


-1.0 


(4.4) 


-6.2 


(-11.4) 


2.1 


(20.4) 


10:1 


0.15 


-16.0 


(-18.4) 


-13.5 


(7.7) 


6.5 


(12.7) 


0.6 


(3.9) 


-18.0 


(-23.8) 


-2.6 


(12.0) 


10:1 


0.40 


-17.5 


(1.4) 


-23.3 


(-7.7) 


5.8 


(11.8) 


-0.2 


(4.4) 


-21.6 


(-7.6) 


-7.2 


(7.1) 



are systems which have remained undisturbed long enough 
for vigorous radiative cooling to have produced a positive 
central temperature gradient and high central gas densities. 

To initialize the structure of our systems, we hav e fol- 
lowed the analytic presc ription of iBabul et. al. I (|2002T ) and 
iMcCarthv et. al] |2004) to produce clusters which conform 
with recent theoretical and observational insights into clus- 
ter structure. The dark matter d ensity profiles of our sys- 
tems follow an NFW-li ke form l|Navarro. Frenk.~fc White! 
I1996J ; I Moore et al.l Il998l ) with the central asymptotic log- 
arithmic slope chosen to be j3 — 1.4 and the scale radius 
(r s ) selected to yield a concentration c = R2oa/r s = 2.6 
(we will use _Ra throughout to indicate the radius within 
which the mean density of the system is A times the critical 
density, p c — 3Hq/8ttG; R200 = 1785 kpc, R 5 oo = 1166 kpc 
and Rcooi = 180 kpc initially for our primary systems). The 
initial gas density and temperature profiles of the clusters 
are set by requiring that (1) the gas be in hydrostatic equi- 
librium within the halo, (2) the ratio of gas mass to dark 
matter mass within the virial radius be f2b/(f2 m — fit), and 
(3) the initial gas entropy 5 scale as S(r) tx r ' over the bulk 
of the cluster body. We normalize the entropy profiles such 
that the temperature of the ICM at R V i r is half the virial 
temperature. 

We construct orbits for our systems which produce spec- 
ified radial and tangential velocities for the secondary sys- 
tem (v r and vt respectively) when its centre of mass reaches 
the virial radius of the primary (R V i r ). For each of the three 
mass ratios we study, we examine three orbits selected to 
produce a typical value of v r (Rvir) and to cover a signif- 
icant range of the transverse velocity Vt(Rvir) giving rise 
to mergers found in cosmological d ark matter simulations 
l|Tormenlll997l ; IVitvitska et al.ll2002l) . 

Throughout our analysis we shall refer to three coordi- 
nate axes; x and y will denote directions in the plane of the 
initial orbit (with the initial separation of the systems be- 
ing along the x axis) and z the direction orthogonal to this 
plane. We shall also distinguish the two systems by calling 



5 We use the standard proxy for entropy given by S = kT/nJ 
with n e and T representing the electron density and temperature 
of the gas. 



the more massive system the "primary" system (we arbi- 
trarily choose one to be the primary in the 1:1 cases, with 
no consequences for our results given their high degree of 
symmetry) and the less massive system as the "secondary" . 

2.2 Dynamical evolution 

In Poole06a we illustrated the generic dynamical progression 
which each of our simulations proceed through, identifying 
five distinct stages: a pre-interaction phase, first core-core in- 
teraction, apocentric passage, secondary core accretion, and 
relaxation. These stages feature prominently in the discus- 
sion which follows so we will briefly review their progression 
here, for those who have not read Poole06a. 

As the cores accelerate towards each other during pre- 
interaction, a pair of shock fronts materialize and are driven 
towards each core, heating and compressing them briefly. 
At t c i ose st the cores of the two systems reach closest ap- 
proach and these effects reach their maximum strength. In 
every case (including the head-on collisions), some part of 
the secondary's core survives its first encounter with the pri- 
mary, tidally stripped into large cool streamers in off-axis 
cases and strings of several clumps in the on-axis cases. At 
tapo the disturbed cores reach maximum separation. After 
subsequently reaching second pericentric passage at t accrete, 
no observable trace of the secondary core remains. This pe- 
riod marks the beginning of a prolonged period of accretion 
when the plume of material generated from the disruption 
of the secondary's core accretes as a high velocity stream. 
This instigates a second episode of core heating followed by 
a period of relaxation. At approximately treiax, the system 
exhibits no obvious substructure in simulated 50ks Chan- 
dra images with the system placed at z = 0.1 (see Poole06a 
for details on how this is done). 

This account of our systems' dynamical evolution is 
terse and we refer readers to Poole06a for a more detailed 
and quantitative account of the various processes involved. 

2.3 Some comments regarding the frequency of 
mergers 

In this paper, we seek to understand the role which merg- 
ers are likely playing in shaping the statistical properties of 
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clusters as studied through observed scaling relations. It is 
important to note however that our simulations certainly 
are not uniformly represented in these statistics. 

In current theories of hierarchical cluster formation, the 
likelihood of merger events is expected to have a strong 
dependence on the mass ratio of the interacting systems. 
Equal mass mergers are expected to be rare while 10:1 
mergers are expected to be common. Recent studies exam- 
ining merger statistics in cos molo gical N-body simula tions 
include ICohn fc White! (|2005h and lWetzel et all (l2006r i. For 
the simulations presented by IWetzel et al.1 (120061 ). the frac- 
tions of clusters which have experienced 2:1, 3:1 and 10:1 
mergers since z — 0.5 are 23%, 49% and 92% respectively 
(private communication; from a sample of 37388 clusters 
with M > 10 14 M Q at z = generated in a ACDM cos- 
mology with Qm = 0.3, h = 0.7, n = 1 and erg = 0.9). 
These statistics are generated using a mas jump criteria 
over conformal time intervals of St = 130/i _1 Mpc. However, 
since multiple mergers which occur over this interval will be 
counted as single events, these numbers are certainly biased 
high. These results are supported qualitatively by an ana- 
lytic extended Press-Schechter analysis we have performed 
on a set of merg e r tree s generated following the approach of 
iTavlor fc Babul (|2004 ). 

Thus, our 3:1 mergers should be representative of the 
most massive merger events statistically important to the 
cluster population. Although there is evidence that equal 
mass merg ers have occurred (e.g. the "Cloverleaf cluster" 
studied by iFinoguenov et all 120051 ) , they are rare and we 
include them in this study primarily as theoretically inter- 
esting limiting cases. The vast majority of clusters should 
have experienced at least one merger similar to our 10:1 
cases. 

For these reasons, when considering the effects of merg- 
ers on scaling relations, we will primarily be concerned about 
the effects of 3:1 and 10:1 events. 



3 LUMINOSITY AND TEMPERATURE 

In this section we present the details of how we compute 
global luminosities and temperatures for our simulations. 
We shall find that the temperatures of our merger remnants 
scale in accord with observed mass scaling relations, but 
that the luminosity does not. Excising the cores of our clus- 
ters improves this discrepancy. We will then illustrate the 
consequences of this for the L x — T x scaling relation. 



3.1 Luminosity 

A cluster merger leads to transient increases in both X- 
ray luminosity (L x ) and temperature (T x ) which can pose 
signif icant problems for cosmological studies (|Randall et al.l 
2002). In their study, RS01 found that after an initial tran- 
sient increase at the time of first pericentric passage, the 
luminosity of the system drops to a fraction of its initial 
value with a second smaller peak ~ 2 Gyrs later when the 
secondary core returns from apocentric passage. These sharp 
transient increases in L x can translate into a skewing of the 
high-mass end of mass functions derived from observed X- 
ray luminosity functions and mass-luminosity scaling rela- 



tions. This can translate into a upward bias in erg of as much 
as 20% and a downward bias in Qm of comparable extent. 

In Fig. [T]we present in red the temporal evolution of the 
bolometric X-ray luminosities of our simulations integrated 
within the varying radius 7 i?5oo(i)- Although detailed com- 
parisons to RS01 are complicated by the fact that we are 
not integrating over the entire simulation volume, we see 
a similar double peaked evolution in luminosity to what is 
noted in their work. Furthermore, we find a similar trend in 
the reduction of the final remnant's luminosity with impact 
parameter. 

This plot illustrates an evolutionary progression with 
features which we will show are common to all of the global 
X-ray and SZ observables we study in this paper. After a 
short period with very little evolution, the system experi- 
ences a dramatic "first transient" increase following t which 
reaches maximum amplitude at t c i osest . This is a product of 
the increased densities and temperatures resulting from the 
compressive forces and shocks formed during the secondary's 
approach to its first pericentric passage through the primary. 
Following t c i 03es t, the system returns to a state similar to it's 
initial conditions. Depending on the observable, it can then 
experience a "second transient" due to the recollapse of the 
system in head-on cases or the second pericentric passage 
of the surviving portion of the secondary core in off-axis 
cases. This second transient is weak for the luminosity of 
our systems but as we will see, can be more significant for 
other observables (e.g. X-ray temperature). Following the 
second transient, the system then evolves quiescently as the 
system relaxes. This progression is in good qualitative agree- 
ment with the results of previous idealized merger studies 
(e.g. RS01). Since all of the X-ray and SZ observables we 
study have direct dependencies on gas density and temper- 
ature, they all evolve with this same qualitative behaviour. 
This leads to a similarly generic qualitative evolution (with 
occasional important differences) in all the scaling relations 
we will examine. 

There is however one notable and generic feature in 
the luminosity evolution of our simulations not present in 
RS01: a persistent increase in luminosity after t re i ax . This 
increase follows the relation L x (t) — L x (t re i ax ) exp (at) 
where a = 0.08 — 0.14 Gyr -1 . We also plot, in blue on 
this figure, the evolution of each system's luminosity with 
the projected central O.I-R200 excised (a common measure 
of the cooling radius of the system; -R200 = 1785 kpcand 
Rcooi = 180 kpc initially for our primary systems). These 
curves (and an examination of the evolving surface bright- 
ness profiles of our simulations) reveal that most of this 
change is arising from increases in the surface brightness 
of the system within this excised region. This is a result of 
increasing central gas densities arising primarily from cool- 
ing and (to a lesser extent) the reaccretion of material to 
the center of the system, dispersed to large radii during the 
interaction. 

Throughout this paper we will examine the degree by 
which our simulated mergers evolve to final states which 
maintain observed scaling relations. We do not necessarily 
expect our simulations to preserve these scaling relations 



7 Throughout this paper, we adopt the peak of the projected 
bolometric X-ray surface brightness as the center of the system. 
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v t /V c =0.00 v/V^O.io v,/V c = 0.40 




t [Gyrs] 

Figure 1. Bolometric X-ray luminosities (L x ) measured within i?5oo(*)- The thin blue curve traces the evolution of the system with 
the central O.I-R200 excised from consideration (i.e. "core corrected" luminosities) while the thick red curve traces the evolution with 
this region included (i.e. "uncorrected" luminosities). Thick horizontal dashes indicate the initial bolometric luminosities scaled by 
the ch ange in total mass within -R500 from initial to final states using the observed Mt — L x scaling relation of Rciprich & Bohringer 
J2002h - Vertical dotted lines indicate (from left to right) the times t (virial crossing), t c i oses t (first pericenter), t aC crete (second 
pericenter) and t re i ax (the moment our remnants would appear relaxed to a 50 ks Chandra exposure at z = 0.1). Black text around 
the boundary indicates the mass ratio and Vt /V c depicted by each panel. 



but we use these measures as a means of estimating the 
degree by which we expect each event to contribute to the 
observed relation's scatter. For this reason, we have plotted 
thick horizontal dashes on Fig. [TJto indicate the bolometric 
luminosities we would expect if our remnant systems evolved 
to states which prese rve the observed Mt — L x s caling re- 
lation determined bv lReiprich fc B ohringer (2002]). We use 
the total mass integrated within R500 for this purpose. In 
head-on cases we see that the total luminosity (in red) of 
our remnant systems generally manage to recover their ex- 
pected values by the end of the simulation. However, the 
luminosities of our off-axis merger remnants are systemati- 
cally lower (by 16 — 50%; see table [TJ than expected from 
the observed mass scaling relation. This discrepancy gener- 
ally increases with impact paramet e r. It is consistant with 
the simulations of iMcCarthv et al.l ((200 7a) who find that 
off-axis mergers produce remnant cores with systematically 
higher entropies and hence, lower densities and luminosities. 



Luminosities computed with the central O.I-R200 excised 
(in blue) are more successful at preserving the observed mass 
scaling relation. The discrepancies of the final luminosities 
from these scalings are only 3 — 28% (see table [TJ in these 
cases. 



3.2 Temperature 

Recently several authors have raised concerns that tradi- 
tional methods of computing ICM temperatures in theo- 
retical studies may be introduc ing systematic biases dur- 
ing comparisons to o bservations (IMathiesen fc Evrardll200ll ; 
iGardini et all |2004| ; iMazzotta et al.ll2004 IVikhlininl 120061 ) ■ 
This bias is introduced because the ICM of both theoreti- 
cal and observed systems is not isothermal and the methods 
used to compute representative temperatures have been fun- 
damentally different. 

Most theoretical measures of temperature are computed 
through one of a variety of possible weighted averages while 
observational temperatures are obtained from fits of plasma 
models to observed spectra. The most widely used theoret- 
ical measure has been the emission-weighted temperature 
computed from 



/ k(T)n 2 e T dV 



(1) 



/A(T)nl dV 

where A(T) is the bolometric emissivity of the ICM and 
n e its electron density. More accurate but similarly conve- 
nient and computationally inexpensive "sp ectroscopic- like" 
weigh tings have recen t ly bee n proposed by Ma zzotta et al.l 
<|2004l ) and IVikhlininl j 2006h. From compari sons to hydro- 
dynamic simulations, IMazzotta et al.l (2004) find that the 
weighting 
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Figure 2. Temperatures (T x ) measured within Rsoo(t) using several techniques discussed in the text. Emission- weighted temperatures 
(Tew) are ill ustrated with red cur ves, spectrally fit temperatures (T spec ) with black curves, and "spectroscopic-like" (T si , computed 
according to Mazzotta et al. 2004) temperatures with blue curves. The thick magenta curves are spectrally fit temperatures with the 
cool cores (r < O.I-R200) excised. Horizont al dashes indica te sp ectral temperatures scaled fr om their initial values using the observed 
Aft — T x scaling relation from the data of iHorner 1 ll200ll) and jRciprich & Bohringcr ( 20oJ . Vertical dotted lines indicate (from left 
to right) the times to (virial crossing), t c i oses t (first pericenter), t aC crete (second pericenter) and t re i ax (the moment our remnants 
would appear relaxed to a 50 ks Chandra exposure at 2 = 0.1). Black text around the boundary indicates the mass ratio and vt/V c 
depicted by each panel. Note the change in vertical scale. 



Jngr 1 — dv 

Jn 2 e T- a dV 



(2) 



with a = 0.75 gives results in much better agreement to tem- 
peratures obtained from simulated Chandra observations. In 
their analysis, they typically find T B i to be lower than T ew 
by 20 — 30%. This relation is calibrated from mixtures of 
isothermal plasmas of two temperatures, with variations of 
T s i from spectral fits being less than 10 %. 

The weighting of IVikhlininl (l2006t ) is somewhat more 
complicated but more accurately takes into account effects 
on T x introduced by instrumental variations and metallicity 
dependent line emission. They fin d a very similar power law 
scaling (a = 0.79 for Chandra) as lMazzotta et al.l (|2004l ) for 
continuum contributions but significant spatial variations 
may exist between the two measures, particularly in cluster 
cores which are dominated by low temperature line-emitting 
gas. However, we are only interested in globally averaged 
temperatures in this work, which are dominated by contin- 
uum emission from hot gas. Under these circumstances, very 
little difference is expected between these measures and so 
w e shall presently f ocus o ur attention on the simpler method 
of iMazzotta et alj (|2004l ). 

Our systems constitute mixtures of gas covering an ex- 
tremely wide range of densities and temperatures, naturally 
raising concerns that the specific weightings determined for 



T s i may breakdown occasionally during the interaction. To 
check this we have also computed spectrally fit tempera- 
tures (T spec ). To produce this measure we compute inte- 
grated spectra over the Chandra bandpass (0.7-7 keV with 
20 eV resolution) and fit isothermal plasmas to them (we 
assume a metallicity Z = 0.3^© for all calculations in this 
paper). This is done via \ 2 minimization weighted by the 
photon flux of each spectral bin. The instrumental response 
of Chandra is included in this analysis as well as Galactic 
absorption at a level of Nh = 2 x 10 20 cm -2 . 

In the analysis of RS01, the emission-weighted temper- 
ature of the system remains relatively constant, interrupted 
only by the pair of transients discussed above, during first 
and second pericentric passages when the cluster cores are 
shock heated and compressed. As their systems relax, they 
recover to slightly higher emission-weighted temperatures in 
equal mass events and to approximately the initial temper- 
ature in 3:1 cases. 

In Fig. [2]we present the temporal evolution of the three 
temperature measures discussed above, integrated within 
Rsooit) for our simulations. In all cases we witness the famil- 
iar double peaked behaviour noted by RS01. We addition- 
ally observe the trend noted by RS01 that the amplitude of 
the initial peak decreases with impact parameter while the 
relative strength of the two peaks becomes more similar. 

It has been noted by IMazzotta et alj l|2004 ) that the 
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Figure 3. X-ray temperatures (spectrally fit) plotted against bolometric luminosities for our simulations (both integrated within R500) 
and compared against observations. Red four and blue five point stars indicate the states of the system at t and t re [ ax respectively. 
The red vector indicates the evolution from t a to t c i oaest while the blue vector illustrates 3Gyr of evolution following t re [ ax . The 
magenta line tracks the evolution from t c i oaeat to t re i ax (generally 4-5 Gyrs) during which the system is visibly disturbed. Green 
cr osses ind i cate t he remnant state necessary to preserve the observed mass-scaling relations. Black points are the observed catalogue 
of iHorner I <200ll) . Dashed lines are our fiducial analytic entropy injection models (So = 10,100,300, and 500keVcm , increasing from 
right to left). Black text around the boundary indicates the mass ratio and vt/V c depicted by each panel. 
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mismatch between T ew and T s i increases with the thermal 
complexity of the system. This is reflected in our simula- 
tions by an increased discrepancy between T £W and our other 
temperature measures during the transient temperature in- 
creases at first and second pericentric passage. This discrep- 
ancy, which is 15 — 40% larger than the difference between 
T em and T s i at all other times, suggests that the transient in- 
creases determined by RS01 (who present emission-weighted 
temperatures) are significantly overestimated. Consequently 
it suggests that the influence of transient te mperature in- 
crease s on measures of as determined by iRandall et al.l 
(2002) are overestimated. 

As expected, T B i and T spsc are systematically lower than 
T ew at all times. Initially, this difference is small (~ 3%) but 
at late times it can be as much as ~ 10%. Given that T s i 
produces results in very good agreement with T apec through- 
out our merger interactions (less than 5% difference in all 
cases and at all times except 500 Myrs during first pericen- 
ter in the 1:1 near-axis cases where we see a 20% difference), 
we conclude that the methods of spectrally weighted tem- 
perature averages developed by iMazzotta et alj (I2004T ) and 
IVikhl ininl (|2006T ) are indeed robust during the complicated 
evolution of a cluster merger. We shall use our spectrally 
fit temperatures for the remainder of our analysis however, 
unless otherwise stated. 

Given that there is such a persistent offset between 
T em and T spec , the normalization of theoretical tempera- 
ture scaling relations generated with this approach are likely 
too high. This has important implications for cosmological 
studies utilizing emission- weighted temperature functions as 
well. 

In Fig. [2] we place thick horizontal dashes to indicate 
the spectral temperature the final remnant would have if 
it were scaled from the initial temperature by the Mt — T x 
relation deriv ed from the temperatures o f Irlorner I (|200ll ) 
and masses o f iReiprich fc Bohringerj (2002). We see that im- 
mediately following t re i ax , the system is significantly cooler 
than this scaled value. Late increases evolve according to 
T x (t) = T x (t relax )exp(at) (with a = 0.01 -O.OSGyr" 1 ), 
taking the system to roughly the observed mass-scaled tem- 
perature by the ends of our simulations in all cases. Of the 
three temperature measures, T em produces results which 
most often and most significantly fail to yield the mass- 
scaled result, suggesting that previous theoretical studies 
of the L x — T x and Mt — T x relations which have employed 
T em as a temperature measure have exaggerated the scatter 
in this relation introduced by mergers. In Table [T] we list 
the discrepancies of our remnants' temperatures from the 
observed mass-scaling relation. 

3.3 L x — T x relation 

In Fig. [3] we present the evolution of our systems in the 
L x — T x plane. We illustrate the evolution from the initial 
state (red 4 point star) to t c i osest (first pericentric passage) 
with a red arrow and the evolution for 3 Gyrs following 
t r ei ax (blue five point star) with a blue arrow. In each case, 
the final L x and T x , scaled fro m their initial values by the ob - 
served mass scaling relations of lReiprich fc Bohringerl (120021 ') 
an orner I (|2001h . are illustrated by a green cross. Magenta 
curves track the evolution of the system from t c i osest to 
trelax, while the system appears significantly disturbed (an 



interval of 4 — 5 Gyrs, except for the 10:1 vt/V c ~ 0.4 case 
which takes 7 Gyrs to appear relaxed) . This format will be 
used for all the scaling relations we subsequently present. 

It has been known for some time that variations in the 
structure of the central regions of clusters contribute signif- 
icantly to the scatter in both luminosity and temperature 
mass scaling relations. For this reason, and out of a desire 
to obtain the tightest relations possible (usually for the pur- 
poses of obtaining mass functions for cosmological studies), 
authors have typically presented "core-corrected" observa- 
tions for studying the L x — T x relation. It has been shown 
by M04 however that there is very interesting structure in 
the L x — T x plane which correlates with the morphology of 
the cores when the central regions are retained for analy- 
sis. To explore the consequences (and possible contributions 
of) mergers to this s tructure, we use the "uncorrected" cat- 
alogue presented by iHorner I (|200lf ) for comparison of our 
results to observations (presented in black on Fig. [3} • 

On Fig. [3] (and most subsequent scaling relation plots) 
we have also placed dashed lines depicting theoretical mass 
scaling relations for a set of fiducial analytic entropy injec- 
tion models for systems with various minimum entropies. 
These models ar e com puted following the procedure of 
iMcCarthv et. qZ~l ((2004) with the following exceptions: we 
use the mass-concentration relation fro m the Mill e nnium 
Simulation, the baseline entropy profile of lVoit et all (|2005h . 
and assume a bias M g (< R2oo)/Mt(< -R200) = O.9f2(,/O m . 
Four lines are plotted for core entropies of 10, 100, 300 and 
500 k eVcm 2 (increasing from right to left). In lBabul et. al. I 
(2002) it was found that the 300keVcm 2 model best fits the 
median relation of the data. This model is thus indicated 
with a thicker line type and we shall consider it to repre- 
sent the median relation of the data in all discussions which 
follow. 

In Section f3. II we discussed the common temporal evo- 
lution of the global X-ray and SZ observables we study in 
this paper. We also noted that this leads to a common qual- 
itative evolution in their scaling relations as well. We see 
the first illustration of this progression in Fig. [3] Initially, 
due to the shocks and compressive forces generated as the 
interacting systems reach first pericenter, the system's lumi- 
nosity and temperature both increase sharply (illustrated by 
red vectors). In all cases, these vectors run roughly parallel 
to the mass scaling relations of our fiducial models in this 
plane, supp orting similar finding s in previous merger stud- 
ies (e.g. iRowlev et. al. I [20041 ). This will be the general 
case (with some important differences) in the other scaling 
relations we shall study as well. 

Following the first transient, the system returns to 
nearly it's initial state at t and subsequently evolves er- 
ratically until t re i ax (when the system appears relaxed in 
our simulated Chandra observations and is virialized within 
-Rsoo)- The magenta curve traces the system's evolution dur- 
ing this period, which lasts 4 — 5 Gyrs. During this time it 
experiences a second transient (illustrated by a sharp jump 
along the fiducial models in the magenta curve) when the 
system recollapses in the head-on cases or the surviving por- 
tion of the secondary's core returns for a second pericentric 
passage in the off-axis cases. The magnitude of this second 
transient increases with impact parameter due to the re- 
duced disturbance of the secondary core. It is during this 
period when most of the displacement across the scatter of 
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Figure 4. Total masses within -R500M computed from our simulations. The actual total mass is plotted as a thick black curve, the 
hydrostatic equilibrium mass (computed from EqnfJ} is plotted as a red curve and the isothermal /3-model mass (computed from 
Eqn|8} is plotted as a dashed blue curve. Vertical dotted lines indicate (from left to right) t (the time of virial crossing), t c i osest (first 
pericenter), taccrete (second pericenter) and t re i ax (the moment our remnants would appear relaxed to a 50 ks Chandra exposure at 
2 = 0.1). Black text around the boundary indicates the mass ratio and vt/V c depicted by each panel. 



the observations occurs. Although difficult to illustrate with 
this plot, much this displacement occurs during the accre- 
tion of streams following the secondary's second pericentric 
passage in several cases. We will examine in mo re detail the 
complic ated processes driving this evolution in lPoole et al.l 
(|2007al ) where we will examine the effects of mergers on the 
cores of our systems. 

Once second pericentric passage has passed and obvi- 
ous substructure dissolved, the blue vector (which illustrates 
3 Gyrs of evolution following t re i ax ) shows that the system 
once again evolves parallel to the mass scaling relations of 
our fiducial models. This occurs as cooling leads to a denser, 
more luminous core and material dispersed by the interac- 
tion reaccretes to the system centre, deepening its gravita- 
tional potential and increasing its temperature. 

The initial conditions we have chosen place our primary 
systems on the high luminosity side of the scatter in the 
L x — T x plane, where systems with compact cool cores are 
typically found (M04). In Section r3.2l we found that mergers 
tend to evolve to states which preserve the observed mass- 
temperature scaling relation but in Section 13.11 found that 
they create remnants which are under luminous with respect 
to the observed mass-luminosity scaling relation. In Fig. [3] 
we see the consequences of these trends: off-axis mergers 
tend to push typical compact cold core systems towards the 
low-luminosity side of the observed L x — T x relation. 

In our off-axis 1:1 cases (which we must emphasize, are 
rarely observed and expected to be exceptionally uncom- 
mon) our mergers are able to drive systems to states which 



cover the majority of the scatter in this plane. However, in 
the more common situations represented by our 3:1 mergers, 
states represented by our fiducial entropy injection models 
with entropies higher than S = 300 keVcm 2 are not reached, 
not even transiently. Our 10:1 mergers fail to reach the me- 
dian relation of the data under any circumstance. 

When this plot is generated utilizing emission-weighted 
temperatures, the full dispersion of the observations is cov- 
ered by both our 1:1 and 3:1 mergers. It is now well es- 
tablished that this temperature measure is not viable how- 
ever, so this does not represent a means by which to ac- 
count for the observed scatter. Instead, it further empha- 
sizes the importance of computing proper spectral temper- 
atures when comparing theoretical models of X-ray clusters 
to observations. Furthermore, not only should the normal- 
ization of past temperature scaling relations generated with 
the emission-weighted measure be viewed with some scep- 
ticism (as argued in Section 1 3 . 2 p , the scatter should be as 
well. 



4 MASS 

To date, all catalogues of cluster masses involving significant 
sample sizes have been compiled from X-ray observations. 
The methods used are thus not direct (such as with lensing 
measures) but are inferred from temperature and surface 
brightness profiles using methods which may be highly sus- 
ceptible to systematic biases, particularly during a merger. 
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These approaches generally assume that the system can be 
described as being in a state of hydrostatic equilibrium; a 
condition given by 



_idP _ GM t (< r) 
Pa dr ~ r 2 

This yields a mass profile of the general form 



(3) 



PA 



(1 + (r/r c ) 2 ) 3/3/2 



(7) 



Under these assumptions, fitting Eqn. [6] to the surface 
brightness profile yields the logarithmic derivative of the 
density profile, producing the following relation for the total 
mass within -R500 



kT(r) 
MM = - v ' r 
Gpm p 



dlogp 9 (r) + dlogT(r) 



dlog r 



dlog 



(4) 



M (R 5OO ) = 



where G and k are Newton's and Boltzmann's constants 
respectively, T(r) is the temperature profile of the system, 
p, the mean molecular weight, m p the mass of the proton and 
p g (r) the gas density profile. In Poole06a we examined the 
evolution of a hydrostatic disequilibrium parameter given by 



H = 1 + 



-1 dp 
Pg dr 



2 GM t (< r) 



(5) 



and illustrated that clusters which appear relaxed can vi- 
olate the condition of hydrostatic equilibrium by H — 
10 — 15% and occasionally even more. How does this affect 
masses obtained from Eqn. [4p 

In Fig. [4] we present (in red) the evolution of the mass 
of our systems within i?soo(i) computed directly from our 
simulations using this relation. For comparison, the actual 
total mass is presented as well (thick black). We can see from 
this plot that despite the limited validity of the hydrostatic 
equilibrium condition, masses computed from Eqn. [5] after 
treiax (when the system appears as a single undisturbed 
remnant) generally reproduce the system mass to within ~ 
5 — 10%, with a slight systematic bias towards higher values. 

To implement Eqn.|4]to measure cluster masses without 
additional assumptions, accurate measurements of p g (r) and 
T(r) must be made. The observational challenges of doing so 
have been considerable but with the increasing availability of 
high quality Chandra and XMM observations, s ignifi cant 
progress has recently been made. IVikhlinin" et all ([20061 ) for 
instance have presented a method of gas density and tem- 
perature deprojection which they have demonstrated to be 
accurate to a few percent for apparently relaxed systems 
jNagai et al]|2006j ). 

Presently however, statistically representative Chan- 
dra and XMM catalogues of sufficient size to study mass 
scaling relations are not available. Fo r this reason, we shall 
use fo r comparison the catalogue of iReiprich fc Bohringerl 
(2002) who report masses within -R500 measured from 
ROSAT and ASCA observations. These authors measure 
masses assuming that their observed systems are in hydro- 
static equilibrium but make the further assumptions that 
the gas is isothermal and that it's density profile follows 
a parameterized form given by the /3-model. For an X-ray 
surface brightness profile given by 



S x {rp) — 



So 



(1 + (r P /r c y)W-^ 2 
the gas density profile of an isothermal gas is given by 



(6) 



3/3kT x R 5 oo (R50o/r c ) 
Gpm p 1 + (R 5 oo/r c ) 2 



(8) 



Temperature and /3-model fi ts in the results presented by 
IReiprich fc Bohringerl {2002) are fit to the flux which ex- 
tends to the observed radius of the system (with the central 
O.I-R200 excised for the temperature fits). The observed ra- 
dius is generally comparable to -R500 which is the radius we 
use to generate the isothermal /3-model masses (computed 
from Eqn. [8} we present in Fig. [4] (in blue). 

From this figure we can see that the isothermal /3-model 
systematically underestimates the mass of our systems by 
25 — 40%. This result is in good agre ement with the find- 
ings of several o t her authors (e.g . lHallman et alj 120061 ; 
iRasia et al] 120061 ; iKav et alj 12004 ). The general trends in 
the scaling relations should not be affected much by this 
bias, but the normalization will be. This could have signif- 
icant implications for cosmological studies conducted with 
isothermal /3-model masses. 



4.1 Mt - T x relation 

In Fig. [5] we present the evolution of our isothermal /3- 
model masses against temperature an d compa r e the m to 
the observed tem perature catalogue of Hornerl (|200lh and 
mass catalogue of IReiprich fc Bohringern 2002T ) . We follow 
the same format here as for Fig. [3j Red four and blue five 
point stars indicate the system's state initially and when 
it appears relaxed at t re i ax . The red vector illustrates the 
transient evolution from t to t c i oses t (first pericentric pas- 
sage) and the blue vector shows 3 Gyrs of evolution fol- 
lowing treiax- The magenta line traces the evolution of the 
system between t c i osest and treiax, when the system looks 
visibly disturbed. Green crosses indicate the /3-model masses 
of our remnants scaled from their initial values by the ac- 
tual change of mass within 7?soo, and their temperatures 
similarly scaled by the observed M t — T x scaling relation. 
Dashed lines indicate our fiducial entropy injection models 
(with the entropy minimum increasing with temperature at 
a constant mass). For this plot, we have reduced the masses 
of these models by 30% to reflect the downward bias in 
isothermal /3-model masses found in Section [4] When we 
do so, we find a good agreement between our fiducial mod- 
els and the observed catalogue. Since the observed catalogue 
and our simulation masses are compiled using the same /3- 
model techniques, they should be directly comparable and 
do not require any such adjustment. 

First, we see that there is a considerably smaller degree 
of scatter (in both the observations and our fiducial models) 
in this plane compared to what is seen in the L x — T x plane. 
This is not a surprise given the strong correlation we expect 
between these quantities: it is the depth of the potential 
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Figure 5. Isothermal /3-model masses plotted against spectrally fit temperatures for our simulations (both integrated within R500) 
and compared against observations. Red four and blue five point stars indicate the states of the system at t and t re [ ax respectively. 
The red vector indicates the evolution from t a to t c i oaest while the blue vector illustrates 3Gyr of evolution following t re i ax . The 
magenta line tracks the evolution from t c i oaeat to t re i ax (generally 4-5 Gyrs) during which the system is visibly disturbed. Green 
cr osses ind i cate t he remnant state nec essary to preserve the observ ed mass-scaling relations. Black points are the observed catalogues 
of iHornerl ll200lh (temperatures) and iReiprich & Bohringcr (2002|) (masses). Dashed lines are our fiducial analytic entropy injection 
models (So = 10,100,300, and 500 ke Vein 2 , increasing from left to right) with the masses reduced by 30% to account for the bias in the 
isothermal /3-model approach. Black text around the boundary indicates the mass ratio and vt /V c depicted by each panel. 
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Figure 6. Isothermal /3-model masses plotted against bolometric luminosities for our simulations (both integrated within -R500) and 
compared against observations. Red four and blue five point stars indicate the states of the system at t and t re [ ax respectively. 
The red vector indicates the evolution from t a to t c i osest while the blue vector illustrates 3Gyr of evolution following t re [ ax . The 
magenta line tracks the evolution from t c [ oaeat to t re i ax (generally 4-5 Gyrs) during which the system is visibly disturbed. Green 
cr osses indicate the remnant st ate necessary to preserve the observed mass-scaling relations. Black points are the observed catalogue 
of iReiprich fc Bohrineerl d2002ll . Dashed lines are our fiducial analytic entropy injection models (So = 10,100,300, and 500kcVcm , 
increasing from right to left) with the masses reduced by 30% to account for the bias in the isothermal /3-model approach. Black text 
around the boundary indicates the mass ratio and vt/V c depicted by each panel. 
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(which is set by the mass) which determines the virial tem- 
perature (and hence, X-ray temperature) of the system. For 
the observations, another (potentially more significant) rea- 
son for reduced scatter is the Mp oc T x dependence present 
in the isothermal /3-model. Any transient increase in T x is 
matched by a proportional increase in Mp, reducing the ob- 
served scatter in this plane. 

Next we see that the qualitative evolution of our merg- 
ers in this plane follows the generic pattern we have dis- 
cussed in Sections 13.11 and 13.31 two transient increases in 
both quantities at t c i osest (first pericenter; illustrated by 
the red vector) and t aC crete (second pericenter; illustrated by 
the magenta curve) which take the system roughly along our 
fiducial mass scaling relations, followed by a period of relax- 
ation (illustrated by the blue vector) during which the sys- 
tem slowly evolves upward along the scaling relation. While 
the system is disturbed (magenta curve), it drifts slightly 
across the scatter in the observations. In both planes, the 
initial transient tends to pull the system towards the low 
temperature side of the scatter (i.e. away from the side 
preferentially occupied by NCC systems). 

There are interesting differences between the evolution 
in the L x —T x and M t ~T x planes however. Most notably, our 
mergers occasionally (in the head-on and Vt/V c = 0.15 1:1 
and 3:1 cases) spend very brief (t < 0.5 Gyrs) periods out- 
side of the bounds of our fiducial models, at the very limits 
of the scatter in the observations. These periods correspond 
to those in Fig. [4] during which we witness strong failures 
of the isothermal /3-model at capturing the total mass of 
the system. We also see a couple of brief (t < 0.5 Gyrs) ex- 
cursions to the high temperature side of the observations 
after t accre t e and before t re i ax in the 1:1 off-axis and 3:1 
vt/V c — 0.15 cases. These are not seen in the L x — T x plane. 
This highlights the vulnerability of isothermal /3-model mass 
estimates to the effects of mergers and the spurious contri- 
butions to scatter in the scaling relations likely being intro- 
duced by its use. 

Such rare differences aside, our simulations suggest that 
mergers have a similar effect on the scatter of the Mt — T x 
plane as on the L x — T x plane, for relatively relaxed systems. 
Our (exceptionally rare) 1:1 off-axis mergers produce relaxed 
states which spread across the majority of the scatter in the 
observations while more common 3:1 and 10:1 mergers do 
not create states which reach the high central entropy side of 
the observations preferentially inhabited by NCC systems. 

Recently. [Baloeh et aD ([2006) have studied the M t -L x 
and M t —T x scaling relations, paying attention to the scatter 
introduced by variations in dark matter halo concentrations 
and formation epochs. They find that, although reasonable 
variations in these properties can account for the observed 
scatter in the M t — T x plane, they can not do so for the 
observations in the M t — L x plane. 

Variations in dark matter halo concentrations are a 
product of system-to-system variations in formation histo- 
ries, which regularly consist of several merger events. The 
fact that our simulations generally remain well within the 
scatter of the Mt — T x plane lends support to the claims of 
iBalogh et ail (|2006l ) that mergers do not produce significant 
scatter beyond their effects in driving variations in forma- 
tion time and dark matter concentration. 



4.2 M t — L x relation 

In Fig. |6]we present the evolution of our isothermal /3-model 
masses against luminosity, following the same format as Fig. 
[5] Once again we see the same qualitative evolution wit- 
nessed in the L x — T x plane: two transients at t c i oaes t (first 
pericentric passage; red vector) and t accre te (second pericen- 
tric passage; magenta curve) followed by a period of relax- 
ation (blue vector) during which the system evolves up the 
mass scaling relations. Furthermore, we see the same brief 
failures of the /3-model mass estimates seen in the M t — T x 
plane, which drive the system briefly (< 0.5 Gyrs) outside 
the region encompassed by our fiducial analytic models. 

The scatter in this plane is much larger and it is clear 
that our simulations fail to produce states of sufficiently low 
luminosity to account for the observations in this plane. 
There is a hint that our systems evolve towards the low- 
luminosity side of the observed scaling relation, but only 
briefly (< 0.5 Gyrs) and during the extreme disruptions gen- 
erated during the head-on cases. In the more common off- 
axis cases, the systems tend to evolve away from this side of 
the scatter when significantly disturbed. 



5 SUNYAEV-ZEL'DOVICH EFFECT 

Quite some time has passed since the original realization 
that the spectral signature of inverse Compton scattering of 
the cosmic microwave background (CMB) by the hot ICM 
of galaxy clusters could be detected as a temperature incre- 
ment in the CMB and used as a pro be of the ICM's structure 
(|Sunvaev fc Zeldovichlll970l . I1972T ). Known as the thermal 
Sunyaev-Zel'dovich (SZ) effect, the observational challenges 
involved in its detection have been significant, making abun- 
dant and reliable data slow to emerge. 

Several sensitive large-area SZ surveys will soon be 
available however and due to the weak redshift depen- 
dence on the mass threshold of a cluste r's SZ detectability 
(~ 10 14 M Q , set by confusion limits; see lHolder et al.ll2007l . 
for more details), catalogues of tens of thousands of clus- 
ters extending to high redshifts will soon be available as an 
incredibly powerful tool for studying cluster evolution and 
cosmology. The potential for mergers to introduce biases in 
these catalogues has a similar, but likely less severe poten- 
tial to introduce s ystematic uncertainties in to the analysis 
of these datasets (|Maiumdar fc Mohrl i2003V Furthermore, 
because the SZ effect is dependent on a line of sight inte- 
gral involving a single factor of density (rather than the p 2 
dependence of X-ray emission studies), SZ observations of 
mergers will be useful as a complementary tool to X-ray 
emission observations for determining the structure of the 
ICM. 

The magnitude of a system's thermal SZ signature as a 
function of position on the sky (denoted by 8) is typically 
quantified by the Compton "y-parameter" given by 

y^) = ^T2 f p ^dl (9) 
m e c z J 

where the integral is of the three dimensional pressure 
Pe(r) = n e (f)kT(f) at position r and is performed along 
the line of sight through the system. 
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Figure 7. The evolution of the central Compton parameter (y ) as a function of time for our simulations. Results computed assuming 
an isothermal /3-models are shown with thin blue curves while direct integrations of our simulations (with densities and temperatures 
within 30kpc held constant) are shown with thick red curves. Horizontal dashes represent the final values expected for the /3-modcl 
results from the mass scalings derived from our fiducial models (see Section 13.11 1. Vertical dotted lines indicate (from left to right) 
t (the time of virial crossing), t c i oaest (first pericenter), t accre te (second pericenter) and t re i ax (the moment our remnants would 
appear relaxed to a 50 ks Chandra exposure at z = 0.1). Black text around the boundary indicates the mass ratio and vt/V c depicted 
by each panel. 
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Figure 8. The evolution of integrated SZ fluxes (S^D^/ /„) computed within -R2500 f° r our simulations as a function of time. Results 
computed assuming an isothermal /3-model are shown with thin blue curves while direct integrations from our simulations are shown 
with thick red curves. Horizontal dashes represent the final values expected from the mass scalings derived from our fiducial models (see 
Section [3. It . Vertical dotted lines indicate (from left to right) t a (the time of virial crossing), t c i oseat (first pericenter), taccrete (second 
pericenter) and t re iax (the moment our remnants would appear relaxed to a 50 ks Chandra exposure at 2 = 0.1). Black text around 
the boundary indicates the mass ratio and vt /V c depicted by each panel. 
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Figure 9. Evolution of the SZ Compton parameter (y ) plotted against temperature for our simulations (both integrated within -R500) 
and compared against observations. Red four and blue five point stars indicate the states of the system at t and t re [ ax respectively. 
The red vector indicates the evolution from t a to t c i osest while the blue vector illustrates 3Gyr of evolution following t re [ ax . The 
magenta line tracks the evolution from t c [ oaeat to t re [ ax (generally 4-5 Gyrs) during which the system is visibly disturbed. Green crosses 
in dicate the remnant state necessar y to preserve the obser ved mass-scaling relations. Black points illustrate the o bserved catalogues 
of iReese et al.l ||2002| . triangles) and lLaRooue et al.l ll2006l . squares) with temperatures taken from [Horner] d200lh . Dashed lines are 
our fiducial analytic entropy injection models (So = 10,100,300, and 500keVcm 2 , increasing from left to right). Black text around the 
boundary indicates the mass ratio and vt /V c depicted by each panel. 
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Figure 10. Evolution of the integrated SZ Compton parameter (Sv, integrated within -R2500) plotted against temperature (integrated 
within -R500) for our simulations and compared against observations. Red four and blue five point stars indicate the states of the system 
at t D and t re [ ax respectively. The red vector indicates the evolution from t a to t c i oaeat while the blue vector illustrates 3 Gyr of 
evolution following t re [ ax . The magenta line tracks the evolution from t c [ oaeat to t re [ ax (generally 4-5 Gyrs) during which the system 
is visibly disturbed. Green crosses in dicate the remnant state necessar y to preserve th e obse rved mass-scaling relations. Black points 
illustrat e the o bserved catalogues of iReese et al.l ||2002| . triangles) and lLaRooue et al. (2006, squares) with temperatures taken from 
iHorner 1 ll200ll) . Dashed lines are our fiducial analytic entropy injection models (So = 10,100,300, and 500keVcm 2 , increasing from left 
to right). Black text around the boundary indicates the mass ratio and vt/V c depicted by each panel. 
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Figure 11. Evolution of the SZ Compton parameter (y ) plotted against bolometric luminosity (integrated within R500) for our 
simulations and compared against observations. Red four and blue five point stars indicate the states of the system at t and t re i ax re- 
spectively. The red vector indicates the evolution from t a to t c i osest while the blue vector illustrates 3Gyr of evolution following 
treiax- The magenta line tracks the evolution from t c ; oses t to t re i ax (generally 4-5 Gyrs) during which the system is visibly disturbed. 
Green crosses indicate the remna nt state necess ary to preserve the obs erved mass-scaling relations. Black point s illustrate the observed 
catalogues of I Reese et al.l l|2002l , triangles) and iLaRoque et al.l j200tl squares) with luminosities taken from iHorner I ll200ll) . Dashed 
lines are our fiducial analytic entropy injection models (S a = 10,100,300, and 500keVcm 2 , increasing from right to left). Black text 
around the boundary indicates the mass ratio and vt/V c depicted by each panel. 
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Figure 12. Evolution of the integrated SZ Compton parameter (Sv, integrated within ifesoo) plotted against bolometric luminosity 
(integrated within -R500) for our simulations and compared against observations. Red four and blue five point stars indicate the states 
of the system at t a and t re i ax respectively. The red vector indicates the evolution from t a to t c ; oses t while the blue vector illustrates 
3 Gyr of evolution following t r£ i ax . The magenta line tracks the evolution from t c ; oses t to t re [ ax (generally 4-5 Gyrs) during which the 
system is visibly disturbed. Green crosses i ndicate the r e mnan t state necessar y to preserve the obse rved mass-scaling relations. Black 
point s illustrate the observed catalogues of iReese et al.l (|2002| . triangles) and lLaRoque et al.l J2006I . squares) with luminosities taken 
from [Horner] ll200lh . Dashed lines are our fiducial analytic entropy injection models (So = 10,100,300, and 500keVcm 2 , increasing 
from right to left). Black text around the boundary indicates the mass ratio and vt/V c depicted by each panel. 



The impact of mergers on relaxed X-ray clusters II 19 



Compton parameters through the centres of systems 
(8 = 0; denoted y ) are commonly reported. However, due 
to the difficulty of correcting observed central decrements 
for the contamination from negative side lobes introduced 
by interferometers, observed values are typically determined 
from fitting an assumed functional form of y(9) (convolved 
with the beam of the instrument) to the observed profile, 
y bs(Q)- The SZ profile commonly assumed is that expected 
for an isothermal /3-model 

v'w - 77-r, y ° (3g-D/ 2 ( 10 ) 

(1 + (r/r c y) y " 

with j3 and r c obtained from X-ray observations of the sys- 
tem. The normalization of the fit of y' (8) to y i, s ( 6) yields the 
obser ved value of y (see lLaRoaue et al.l 120061 : iReese et al.l 
2002, for more detailed accounts of this procedure). 

This approach will tend to capture the system's struc- 
ture outside of r c and will treat the core as uniform. Thus, 
since nearly all systems have r c > 30 kpc, variations within 
30kpc due to complex AGN processes will tend to be di- 
minished. In Fig.[7]we demonstrate this fact by plotting the 
time evolution of y (measured in the z-projection) for each 
of our simulations. The thin blue curve illustrates y com- 
puted using the /3-model approach presented above, and the 
thick red curve presents a direct integration of our simula- 
tions with the ICM's central (r < 30 kpc) temperature and 
density set to it's value at 30 kpc. The two results are in 
very good agreement: the /3-model approach often used in 
observational studies does a very good job of capturing the 
structure of the system outside of 30kpc over the entire 
duration of our mergers. 

We estimate that the contribution to y from mate- 
rial within 30 kpc is ~ 20%. Hence, large fluctuations due 
to AGN activity in this region can thus have a significant 
effect on y . Thus, the tendency of the /3-model approach 
to diminish variations in y has both positive and negative 
consequences. On the one hand, reducing the sensitivity of 
y a to regions whose structure is not well understood and 
whose evolution is affected by stochastic effects generated 
by central AGN should lead to much tighter scaling rela- 
tions. This is useful for cosmological studies where measures 
which accurately characterize the total mass of the system 
are desired. On the other hand, reducing the sensitivity of 
this measure to the cluster's most central regions reduces its 
utility for studying the uncertain structure of cluster cores 
and the processes active within them. 

We have also placed thick dashes on Fig. [7] to repre- 
sent the expected mass scaling of y from its initial values 
(we assume 3/0,2/3/0,1 = (AI2/Mi) a , with a determined to 
be 0.93 from our fiducial models since adequate empirical 
measures are not yet available). This is performed for the 
/3-model results which include the central core. We can see 
that the expected mass scaling of y is recovered in each 
case to within 30%, with a systematic tendancy for y to 
fall below this val u e with increasing impact parameter. 

iBenson et aD Hflflj ) have pointed out that y can be 
significantly biased by uncertainties in the assumed spatial 
model of the ICM. They present an analysis of integrated 
Compton fluxes which they find to be more robust to uncer- 
tainties in the assumed form of y'(d). Integrated SZ fluxes 
are computed from Eqn. [9] according to 



y in t{<0) = 2ir f V (O')0'<W (II) 
Jo 

However, SZ measurements are usually expressed in terms 
of frequency dependent flux densities which are given by 



S v = j v (x)y int (12) 

where x = hu/kT is a dimensionless frequency and j v is a 
function which describes the spectral shape of the SZ dis- 
tortion. This function is given by 



>0) = (fee) 2 ^ ' ^ ' 

where 



To facilitate comparisons to observations, which are 
performed at many different frequencies and for systems 
at various redshifts, we remove the frequency and distance 
(redshift) dependencies from our computed flux densities by 
reporting the quantity S V D\/ f v . 

In Fig. [8] we present the evolution of our merging sys- 
tems' SZ flux integrated within a projected radius of -R2500 
(measured in the z-projection). Once again we illustrate two 
quantities in each case: the integrated SZ flux computed di- 
rectly from the simulations and from integrating the best 
fit of y'(6) to y oos (9). We see that integrating the best-fit 
/3-model does an excellent job of tracking the actual value. 
In addition, we have placed thick dashed lines to indicate 
the expected mass scalings of S„ from its initial values (we 
assume S v ^/S v ^ = (M2/Mi) a , with a determined from our 
fiducial models to be 1.63, once again, because adequate em- 
pirical measures are not yet available). We see that this SZ 
measure recovers the mass scaled values far better than y . 

In all our simulations the evolution of both S v and y 
follows the familiar pattern seen for L x (t), T x (t) and M t (t): 
a rapid rise that begins shortly before t c i osest , a subsequent 
drop to values below the initial state, a second (but sub- 
stantially weaker) peak at t accre t e , & slow decrease until 
treiax and then a final rise continuing to the end of the 
simulation. Radiative and adiabatic cooling (arising from 
an expansion of the core) results in a significant reduction 
from the initial value following t accrete- As active cooling 
reestablishes in the remnant cores, the system's integrated 
SZ signature begins increasing. 

For the remainder of this section, we shall examine the 
scaling relations constructed from these quantities and the 
X-ray temperature and luminosity. Observations for com- 
parison to our simulations are still very sparse and statisti- 
cally incomplete at the masses probed by our simulations, 
and we present these plots for completeness and to establish 
theoretical expectations for understanding larger data sets, 
once they become available. 
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5.1 SZ — T x scaling relations 



5.2 SZ— L x scaling relations 



The scaling behaviours of y and S v with T x are presented 
in Figs. l9l and [101 following the format of our previous scal- 
ing relation plots: red four and blue five point stars mark 
the system's state at t D (virial crossing) and t re i ax (when 
the system appears relaxed and is virialized) respectively, 
while red vectors illustrate the system's evolution from t 
to t c i osest (first pericentric passage) and blue vectors mark 
3 Gyrs of evolution following t re i ax . The magenta curve il- 
lustrates the evolution of the system while it is visibly dis- 
turbed between t c i osest and t re i ax . During this time, the sys- 
tem experiences a second transient roughly along the scaling 
relation and drifts slightly across its scatter. 

For observational comparison in our SZ analysis, we use 
th e heterogeneous ca talogues of massive sy stems presented 
by Reese et"afl l|2002l . black triangles) and lLaRoque et al.l 
(2006, black squares). The observed values of S v we present 
are produced from the values of y and the /3-model param- 
eters tabulated by these authors since values of S v are not 
presented in either catalogue. Given the excellent recovery 
of Sv from integration of Eqn. [10] illustrated above, these 
values should be an accurate representation of the actual 
values for these systems. 

Examining the y — T x plane first, we note two general 
trends in our simulations: when the system is most disturbed 
at t c i osest , it tends to be deflected towards high values of y 
and once relaxed, it evolves very little. Since y oc p g T x , the 
first of these observations can be understood as a product 
of the additional increases in y due to higher gas densities 
generated by compression. The second can be understood as 
a consequence of the /3-model treating all systems as con- 
stant density isothermal cores within r c . Even though the 
system may be relaxing and central pressures increasing as 
cooling recovers following t re i ax , the isothermal /3-model will 
not catch this behaviour and very little evolution is seen as 
a result. 



In Figs.[TT]and[T2]we present the evolution of our simulations 
in the y — L x and S v — L x planes, following the same format 
as the SZ— T x scaling relations presented in Section al! 

Examining the y — L x plane first, we immediately no- 
tice the remarkably small variation in the fi ducial mass scal- 
ing rel ations with S , first predicted by iMcCarthv et al.l 
l|2003al ). Our fiducial models capture the median trend of 
the observations and we see that in all cases, our simulated 
systems begin on this predicted scaling relation (red four 
point star) and by t re i ax (blue five point star) have returned 
to it and are evolving along it. 

During first pericentric passage (red vector) and while 
significantly disturbed between t c i oae st and t re i ax (magenta 
curve), the system is over luminous relative to our fiducial 
models. At no time do our mergers evolve to states which are 
significantly under luminous relative to the median trend of 
the observations. However, although the scatter in the ob- 
servations is small, many systems deviate from the fiducial 
relation by a statistically significant amount. Those systems 
which are significantly over luminous relative to our fidu- 
cial models can potentially be accounted for through merger 
activity, but those which are under luminous can not. In- 
creased pressure from AGN activity could be responsible for 
the deviant state of such systems. 

In Fig. [12] we find that the degeneracy in the SZ—L X 
mass scaling relations with S is broken when S„ is used, 
with both our fiducial models and the observations ex- 
hibiting significantly more scatter. Our simulations show 
that while disturbed, the system can swing wildly to states 
both over luminous with respect to the observations and 
our fiducial models (during first pericentric passage) and 
briefly to the lowest luminosity side of the scatter. Once 
relaxed however, states represented by our fiducial models 
with So > 300 keVcm 2 are not produced. 



The range of central entropies chosen for our fiducial 
model were selected to encompass the majority of observed 
systems in the L x — T x plane and they generally do so in the 
other scaling relations we have studied so far. However, in 
the y — T x plane (and all the following SZ scaling relations 
we present) there is a slight bias in the observations towards 
our high-S'o fiducial mod els. This trend has bee n noted in 
previous studies as well ^McCarthy et al.l l2003bh . Without 
knowing the degree to which the observed catalogue is statis- 
tically representative, it is difficult to comment any further 
on this. 

In Fig. [12] we present the evolution of our simulations 
in the S„ — T x plane. Our fiducial models do a better job 
of embracing the scatter in the observed catalogues in this 
case, but there is still a slight bias towards the higher-So 
models. Both the observed scatter and the scatter expected 
from our fiducial models is considerably less. 

In both planes, we once again see that mergers can dis- 
place a system to approximately the So = 300 keVcm 2 fidu- 
cial track, but not to the fiducial tracks of higher-S models 
where many observed systems are located. 



6 THE X-RAY SUNYAEV-ZEL'DOVICH 
PROXY 

Lastly, we have examined the evolution of the X-ray proxy 
for SZ fluxes introduced by iKravtsov et al.l (|2006n . This 
proxy is given by 



Y x (<R)=M g (< R)T X 



(15) 



where M g is the gas mass within R and T x is the X-ray 
temperature measured within R, with the central O.I-R200 
excised. In Fig.[T3]we present the evolution of Y x integrated 
within i?5oo as measured directly from our simulations (red) 
and our isothermal /3-model fits (blue; in the z-projection 
in all cases). The evolution of this quantity is (as expected) 
very similar to that of the integrated SZ flux with a similar 
underestimate of ~ 15 — 25% by the isothermal /3-model 
from actual values. We have once again placed thick dashes 
to indicate the expected mass scalings (derived from our 
fiducial models) from initial values. 
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Figure 13. The X-ray SZ proxy (Y x = M g T x ) for our simulations as a functi on of time. Thick red c urves illustrate values integrated 
directly from our simulations while thin blue curves are computed according to iKravtsov et al. Horizontal dashes represent the 

final values expected from the mass scalings derived from our fiducial models (see Section l3.1l l. Vertical dotted lines indicate (from left 
to right) t (the time of virial crossing), t c i oseat (first pericenter), t accre te (second pericenter) and t re [ ax (the moment our remnants 
would appear relaxed to a 50 ks Chandra exposure at z = 0.1). Black text around the boundary indicates the mass ratio and vt/V c 
depicted by each panel. 




Figure 14. Evolution of the X-ray SZ proxy (Y x ) introduced by IKravtsov et aL (2006) plotted against total mass for our simulations 
(both integrated within -Rsoo)- Red four and blue five point stars indicate the states of the system at t a and t relax respectively. The red 
vector indicates the evolution from to to t c i oaeat while the blue vector illustrates 3 Gyr of evolution following t re i ax . The magenta line 
tracks the evolution from t c i oaf , st to t re i ax (generally 4-5 Gyrs) during which the system is visibly disturbed. Green crosses indicate 
the remnant state necessary to preserve the observed mass-scaling relations. Black text around the boundary indicates the mass ratio 
and vt /V c depicted by each panel. 
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6.1 M t - Y x relation 

In Fig. [14] we present the scaling of the X-ray SZ proxy 
with total actual projected mass within -Rsoo- Astonish- 
ingly, we find that our system remains well confined within 
the 8% scatter bound s of the mean relation determined by 
iKravtsov et al.l l|200fj| ) at all times. There is a slight ten- 
dency for the final remnant to evolve on the low-Kc side of 
the relation while visibly disturbed (magenta line) and to 
evolve to states shifted slightly to the high-F^ side of the 
mean following t re i ax (blue five point star). 

We present results for only the z projection in Fig. [T3] 
but have confirmed that the system remains within the 8% 
scatter bounds at all times in all projections. 



7 DISCUSSION 

Scaling relations are powerful tools for constructing mass 
functions for cosmol ogical studies. However, it has bee n 
known for some time (jFabian et al.ll 19941 ; lMarkevitchlll998l ). 
and we have illustrated in this paper, that variations in core 
properties contribute significantly to the scatter in both X- 
ray and SZ properties. This is undesirable for cosmological 
studies (motivating many authors to exclude the central re- 
gions of clusters from their analysis) but affords an opportu- 
nity to understand the processes which shape cluster cores. 
Such processes include, but are not restricted to, the effects 
of mergers, AGN and various other processes (such as pre- 
heating) which may have contributed to modifying the cen- 
tral entropy (and hence, morphology) of cluster cores. It is 
important to keep this dual role of cluster scaling relations in 
mind when considering global cluster properties: both core 
corrected and uncorrected values are of great utility. 

Keeping the cores in our analysis, we have examined 
scaling relations between luminosity, temperature, mass and 
SZ properties. In most cases, we find that relaxed states 
which cover the full dispersion in the scaling relations can 
only be achieved through extremely rare equal-mass, off-axis 
mergers. Somewhat more common 3:1 mergers are rarely 
able to generate systems (neither transiently nor as relaxed 
remnants) which account for the dispersion in the observed 
scaling relations. In order for mergers to move the system 
the rest of the way across the relation, an additional merger 
of comparable extent would thus be required. To have a sig- 
nificant effect on the statistics of the population this would 
have to occur regularly, requiring the compounding of in- 
frequent events. In short, mergers between typical relaxed 
compact cool core systems can not account for the disper- 
sion in observed scaling relations. 

Ultimately, detailed theoretical studies of the scatter 
in scaling relations must be conduc ted with cosmolo gical 
hydrodynami c simu lations. To date, iKav et all l|2006l) and 
lO'Hara et al.l l|2006l) are the only two cosmological N-body 
studies which have specifically examined the scatter in scal- 
ing relations. Both of these studi es support the results of 
our analysis: lO'Hara et al.l ([2006) find that their simula- 
tions clearly fail to produce enoug h scatter to account for 
observations and IKav et all (|2006h find that nearly all of 
their clusters, regular or irregular, have cool cores. Since ap- 
proximately half the scatter in observed sc aling relations is 
due to NCC systems with warm cores fsee iMcCarthv et al.l 



2007b), we infer that their simulations fail to reproduce the 
observed scatter as well. Furthermore, it is also worth noting 
that the study we have presented in this paper will be useful 
as a tool for interpreting the effects of individual events for 
future scaling relation studies conducted with cosmological 
simulations. 

We have initialized our systems to represent common 
compact cool core systems observed by Chandra and to 
represent a significant range of orbital properties found in 
cosmological dark matter simulations. As a consequence, our 
systems should represent a significant range of typical clus- 
ter merger events. A natural question to ask however, is 
whether systems with initially higher central entropies could 
be driven by mergers to states which account for the disper- 
sion in the scaling relations. 

Such a scenario is effectively a multiple heating event 
scenario in which a system's central entropy becomes el- 
evated (either by preheating, a previous merger event or 
post-collapse AGN activity) and then additionally heated, 
before it is able to relax, to states which single mergers 
between relaxed systems can not generate. This is appeal- 
ing on the grounds that mergers involving heated cores will 
generate lower densities and stimulate less radiative cool- 
ing, tipping the balance between heating and cooling to- 
wards heating. Furthermore, previously heated cores will 
have low er densities, incre asing the efficiency of entropy pro- 
duction (|Voit et al.ll2003t ) from shocks. 

One appealing mechanism (which we shall examine in a 
future study) of generating frequent multiple heating events 
could be the stimulation of AGN activity through mergers. 
There are several stages during the evolution of our mergers 
when gaseous material is driven towards the cores of our 
systems. These times generally occur during the interactions 
of the merging cores (at t c i osest or t accre te) or during the 
accretion of the streams generated from the disruption of 
the secondary's core. This may provide a natural means of 
generating multiple heating events in close succession (and 
thus reproducing the full dispersion of the scaling relations) 
within the frame work of hierarchical clustering. 

Alternatively. iBurns et all (120061 ) have proposed a plau- 
sible means by which cluster mergers (albeit between un- 
relaxed systems at high redshift) could account for the ob- 
served dispersion in cluster properties. Using a set of cosmo- 
logical adaptive mesh refinement (AMR) simulations, they 
argue that early major mergers can prevent the establish- 
ment of cold cores, leading to the production of systems with 
high central entropies. Although this result is preliminary 
and has not yet been reproduced by other authors (partic- 
ularly those utilizing SPH algorithms), it is an interesting 
possibility worth considering. 

Impending observational progress should contribute sig- 
nificantly towards resolving these issues. We have shown 
for example that using the isothermal /3-model to measure 
masses not only leads to a downward bias of ~ 30%, but 
likely also contributes significantly to the observed scatter 
in mass scaling relations. However, it has been demonstrated 
(so far only on limited samples; e. g. IVikhlinin et al.l I2006T ) 
that significant improvements in X-ray mass estimates are 
feasible with high quality X-ray observations using modern 
reduction techniques. Once applied to a large representative 
sample of clusters, interpretation of observed mass-scaling 
relations will be greatly simplified. 
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Furthermore, improvements in available SZ catalogues 
will be extremely valuable, particu larly when combined w ith 
X-ray and weak lensing maps (see iMahdavi et"ai"1l2006al . for 
a discussion of the benefits such a cross correlation will af- 
ford). 



8 SUMMARY 

Building on the result s of previous authors (most notably, 
iRicker fc Sarazi nl l200ll ). we have identified a generic evolu- 
tionary sequence for X-ray and SZ observables during cluster 
mergers. In all the cases we have studied (which we believe 
represent a significant range of cluster merger initial con- 
ditions), these properties (specifically L x , T x , M t , y and 
S v ) exhibit large increases during core-core interactions at 
tciosest and t accre te with an intervening drop from ini- 
tial values due to the expansion of the primary system in 
response to the initial core-core interaction. Our systems 
then experience a subsequent rise in all quantities following 
treiax due to the reestablishment of efficient radiative cool- 
ing in the remnent cores and the reaccretion of displaced 
material. This increase lasts until the conclusion of our sim- 
ulations. 

Because of this common evolution in the global prop- 
erties we have studied, merger systems tend to evolve in a 
generic way in scaling relations constructed from them: there 
is an initial transient which evolves roughly along the mass 
scaling relations, followed by a drift across the relation until 
treiax, after which the system quiescently evolves along the 
mass scaling relation. 

Our study indicates that typical cluster mergers fail to 
generate states which can account for the scatter of observed 
scaling relations. They can push relaxed clusters to the me- 
dian position of the observed scaling relations, but not past 
it (not even transiently). 

Several other, more specific findings of our study, in- 
clude: 

• The luminosity integrated within i?soo(i) recovers the 
observed mass scaling relations only in the head-on cases. 
As mergers are moved off-axis, the remnant's luminosity be- 
comes steadily less than expected from the observed mass 
scaling (discrepancies are 16 — 50%). This discrepancy is still 
present but is reduced to levels of 3 — 28% if the central re- 
gions (r < O.I-R200) are excised. 

• We have studied the evolving temperatures of our merg- 
ers using 3 temperature measures: emission weighted tem- 
pe ratures (T ew ), the "spe ctroscopic- like" average presented 
by iMazzotta et ail (|2004r ) (T 3 i) and isothermal fits to inte- 
grated spectra (T 3pec ). We find that T ew is higher than the 
other two measures at all times and that T 3pec produces re- 
sults which agree with T 3 i to within 5% at nearly all times. 
We thus confirm the utility of this method for generating 
globally averaged temperatures from massive simulated clus- 
ters. 

• T ew yields remnant temperatures significantly in excess 
of what is expected from the observed mass scaling relation 
while the other measures recover the observed mass scaling 
result in all cases to within 15%. This should lead to signifi- 
cantly less scatter in theoretical L x — T x relations computed 
with T 3 i or T 3pec . 



• The temperature spikes following first and second peri- 
centric passage (at t c i 03e3t and t accre t e respectively) are 
15 — 40% larger for the T em measure than the other two 
methods. This suggests that the estimates of the effect 
of merger-driven temperaure spikes on mass functions de- 
termi ned from temperatu re scaling relations utilizing T ew 
(e.g. iRandall et al.ll2002l ) have been exagerated. 

• As a consequence of the mass scaling behaviours of L x 
and T x , head-on mergers ultimately move systems along the 
observed L x — T x scaling relation while off-axis mergers scat- 
ter systems towards lower luminosities. 

• Isothermal /3-models track the evolution of the actual 
mass of the system before the interaction and while appear- 
ing relaxed afterwards, but with a systematic underestimate 
of 25 — 40% supporting the findings of several previous stud- 
ies. 

• The isothermal /3-model does a very good job of re- 
producing the actual results for the SZ observables y and 
S v (< i?250o)- Theoretical mass scaling relations are pre- 
served much better by the integrated measure. 

• R emarkably little scatte r is expected in the y — L x 
plane (|McCarthv et al.ll2003al fbh. Mergers push systems only 
to the high luminosity side of the scatter and the system re- 
laxes back to the median relation afterwards. These results 
make it difficult to use mergers to account for under lumi- 
nous systems in this plane. 

• Regardless of p rojection effects, the SZ proxy of 
iKravtsov et al.l (|2006l ) obeys their reported mass scaling to 
within 8% over the entire duration of all our merger events. 

We conclude that the full dispersion in scaling relations 
can not be produced from simple two body mergers between 
common relaxed compact cool core systems. However, some 
displacement across each plane does occur and this is driven 
primarily by changes in the structure of a cluster's core. 
Thus, the outer (i.e. r > O.I-R200) structure of clusters are 
preserved following merger events. This is satisfying in light 
of recent observations suggesting that the outer structure of 
clusters exhibit a universal structure. 

Mergers certainly alter the structure of compact cool 
cores (although not sufficiently to account for the observed 
range of system properties) however, and we shall explore 
the details of these effects and the processes driving them in 
the next paper of this series. 
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